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• Abstract. Inspiralling compact binaries are ideally suited for application of a high-order post-Newtonian (PN) gravitational 

' wave generation formalism. To be observed by the LIGO and VIRGO detectors, these very relativistic systems (with orbital 

I velocities v ~ 0.5c in the last rotations) require high-accuracy templates predicted by general relativity theory. Recent 

, calculations of the motion and gravitational radiation of compact binaries at the 3PN approximation using the Hadamard 

self-field regularization have left undetermined a few dimensionless coefficients called ambiguity parameters. In this article 
we review the application of dimensional self-field regularization, within Einstein's classical general relativity formulated in D 
space-time dimensions, which finally succeeded in clearing up the problem, by uniquely fixing the values of all the ambiguity 
parameters. 

^ , 

> 1, INTRODUCTION 

oo 

' The problem of the motion and gravitational radiation of compact objects in post-Newtonian (PN) approximations 

\^ . of general relativity is of crucial importance, for at least three reasons. First, the motion of objects at the IPN 

' level,' according to the Einstein-Infeld-Hoffmann equations |1], is routinely taken into account to describe the 

\^ , Solar System dynamics. Second, the gravitational radiation-reaction force, in reaction to the emission of gravitational 

■ radiation, which appears in the equations of motion at the 2.5PN order ^ 1 /c^, has been experimentally verifiedjr^ the 

Q I observation of the secular acceleration of the orbital motion of the Hulse-Taylor binary pulsar PSR 1913H-16 |2, 3, 4, 5]. 

O^' Last but not least, the forthcoming detection and analysis of the gravitational waves emitted by inspiralling compact 

5^ I binaries — two neutron stars or black holes driven into coalescence by emission of gravitational radiation — will 

bX). necessitate the prior knowledge of the equations of motion and radiation field up to high post-Newtonian order 

Inspiralling compact binaries are extremely promising sources of gravitational waves for the detectors LIGO, VIRGO, 
GEO and TAMA. The two compact objects steadily lose by gravitational radiation their orbital binding energy; as 
a result, the orbital separation between them decreases, and the orbital frequency increases. The frequency of the 
gravitational-wave signal, which equals twice the orbital frequency for the dominant harmonics, "chirps" in time {i.e. 
the signal becomes higher and higher pitched) until the two objects collide and merge. 

Strategies to detect and analyze the very weak signals from compact binary inspiral involve matched filtering of a set 
of accurate theoretical template waveforms against the output of the detectors. Several analyses 1 6.l7llMl^ flflfllll have 
shown that, in order to get sufficiently accurate theoretical templates, one must include post-Newtonian effects iroto the 
3PN level at least. To date, the templates have been completed through 3.5PN order for the phase evolution 1 12. flSi. 14,1 . 
and 2.5PN order for the amplitude corrections Spin effects are known for the dominant relativistic spin-orbit 

coupling term at L5PN order and the spin-spin coupling term at 2PN order and also for the 

next-to-leading spin-orbit coupling at 2.5PN order |23, 24, 25, 26]. 

The main point about modelling the inspiralling compact binary is that a model made of two structureless point 
particles, characterized solely by two mass parameters nii and 1112 (and possibly two spins), is sufficient. Indeed, most 
of the non-gravitational effects usually plaguing the dynamics of binary star systems, such as the effects of a magnetic 
field, of an interstellar medium, and so on, are dominated by gravitational effects. However, the real justification for a 
model of point particles is that the effects due to the finite size of the compact bodies are small. Consider for instance 
the influence of the Newtonian quadrupole moments Qi and Q2 induced by tidal interaction between two neutron stars. 



' As usual hPN refers to the terms of order (v/c)^" where v is the typical orbital velocity of the objects and c is the speed of Ught. We shall generally 
abbreviate (v/c)^" as 1/c^". 



We suppose that the neutron stars have no intrinsic spins. Let ai and 02 be the radius of the stars, and L the distance 
between the two centers of mass. We have, for tidal moments. 



where ki and k2 are the star's dimensionless (second) Love numbers, which depend on their internal structure, and are, 
typically, of the order unity. On the other hand, for compact objects, we can introduce their "compactness", defined by 
the dimensionless ratios 

Gnii Gm2 
Ki = T , K2 = 2 , (2) 

which equal ~ 0.2 for neutron stars (depending on their equation of state). The quadrupoles Qi and Q2 will affect both 
the Newtonian binding energy of the two bodies, and the emitted total gravitational wave flux as computed, say, 
using the standard Einstein quadrupole formula. It is known that for inspiralling compact binaries the neutron stars 
are not co-rotating because the tidal synchronization time is much larger than the time left till the coalescence. The 
best models for the fluid motion inside the two neutron stars are the so-called Roche-Riemann ellipsoids |27], which 
have tidally locked figures (the quadrupole moments face each other at any instant during the inspiral), but for which 
the fluid motion has zero circulation in the inertial frame. In the Newtonian approximation we find that within such a 
model (in the case of two identical neutron stars) the orbital phase, ^ = f (Odt where (O is the orbital frequency, reads 

'^'-'^'o=^{i+™KI)'}' 

where x = (Gmo/c^)^/^ is a standard dimensionless PN parameter ^ and where k is the Love number and 
K is the compactness of the neutron star. The first term in the right-hand-side (RHS) of Q corresponds to the 
gravitational-wave damping of two point masses; the second term describes the finite-size effect, which appears as 
a relative correction, proportional to (x/K)^, to the latter radiation damping effect. Because the finite-size effect is 
purely Newtonian, its relative correction ^ {x/K)^ should not depend on c; and indeed the factors cancel out 
in the ratio x/K. However, the compactness K of compact objects is by Eq. of the order unity (or, say, a few 
tenths), therefore the it contains should not be taken into account in order to find the magnitude of the effect 
in this case, and so the real order of magnitude of the relative contribution of the finite-size effect in Eq. (|3} should 
be given by x^ alone. This means that for non-spinning compact objects the finite-size effect should be comparable, 
numerically, to a post-Newtonian correction of magnitude x^ ~ namely 5PN order. ^ This is a much higher 

post-Newtonian order than the one at which we shall investigate the gravitational effects on the phasing formula. 
Using k' = (const A;) ~ 1 and K ~ 0.2 for neutron stars (and the bandwidth of a VIRGO detector between 10 Hz and 
1000 Hz), we find that the cumulative phase error due to the finite-size effect amounts to less that one orbital rotation 
over a total of ^ 16000 produced by the gravitational-wave damping of point masses. The conclusion is that the 
finite-size effect can in general be neglected in comparison with purely gravitational-wave damping effects. ^ Thus the 
appropriate theoretical description of inspiralling compact binaries is by two point masses within the post-Newtonian 
approximation. 

Our strategy to obtain the motion and radiation of a system of two point-like particles at the 3PN order is to start with 
a general form of the 3PN metric, that is valid for a general continuous (smooth) matter distribution. Applying such 
metric to a system of point particles, we find that most of the integrals become divergent at the location of the particles, 
i.e. when x ^ yi (t) or y2{t), where yi (f ) and y2(f ) denote the two trajectories. Consequently, we must supplement the 
calculation by a prescription for how to remove the "infinite part" of these integrals. At this stage different choices for 
a "self-field" regularization, which will take care of the infinite self-field of point particles, are possible. Among them: 

1 . Hadamard's self-field regularization, which has proved to be very convenient for doing practical computations (in 
particular, by computer), but suffers from the important drawback of yielding some ambiguity parameters, which 
cannot be determined within this regularization, at the 3PN order; 



- This result can be derived in the context of relativistic equations of motion, and yields a proof of the so-called "effacement" principle in general 
relativity, according to which the internal structure of the compact bodies does not show up in their motion and emitted radiation which depend only 
on the masses |28]. 

^ But note that for non-compact or moderately compact objects (such as white dwarfs) the Newtonian tidal interaction dominates over the radiation 
damping. 



2. Dimensional self-field regularization, an extremely powerful regularization, that is free of any ambiguities (at 
least up to the 3PN level), and as we shall see permits to uniquely fix the values of the ambiguity parameters 
coming from Hadamard's regularization. However, dimensional regularization will be implemented in the present 
problem not in the general case of an arbitrary space-time dimension D but only in the limit where D ^ 4. 

Dimensional regularization was invented by 't Hooft and Veltman as a mean to preserve the gauge 

symmetry of perturbative quantum field theories. Our basic problem here is to respect the gauge symmetry associated 
with the diffeomorphism invariance of the classical general relativistic description of interacting point masses. Hence, 
we use dimensional regularization not merely as a trick to compute some particular integrals which would otherwise 
be divergent, but as a powerful tool for solving in a consistent way the Einstein field equations with singular point- 
mass sources, while preserving its crucial symmetries. In doing this, we implicitly assume that the correct theory is 
the Einstein general relativity in D space-time dimensions (Section|5}- 

Earlier work on the equations of motion of point masses at the 2PN approximation level fl^ was based on the 
Riesz analytical continuation method | 33J, which consists of replacing the delta-function stress-energy tensor of point 
particles by an auxiliary, smoother source defined from the Riesz kernel, depending on a complex number A (in a 
normal 4-dimensional space-time). However, it was mentioned ^28^1 that the same final result would be obtained by 
considering ordinary delta-function sources but in a space-time of complex dimension D — 4- — A. In other words, 
the Riesz analytic continuation method is equivalent to dimensional regularization. It was also noticed ll28ll that the 
generalization of the Riesz continuation method to higher post-Newtonian orders is not straightforward because of the 
appearance of poles, proportional to A^' — (4 — D)^^ at the 3 PN order 

In the meantime all calculations were performed using the more rudimentary Hadamard regularization |34, 35, 36l 
IstIi . yielding almost complete results at the 3PN order, i.e. complete but for a few ambiguity parameters, which turn 
out to be in fact associated with the latter poles at the 3PN order. Then it was shown |38| how to use dimensional 
regularization within the ADM canonical formalism of general relativity at the 3PN order for the problem of the 
equations of motion. Further work using dimensional regularization |39, 14. .40.1 , finally determined the values of all 
the ambiguity parameters, both in the 3PN equations of motion (parameter called A) and in the 3PN gravitational 
radiation field (parameters ^, K and Q. Further in the same context, we quote Ref. [41] for an alternative approach, 
based on a Feynman diagram expansion, showing how to renormalize using dimensional regularization. 



2. EINSTEIN'S FIELD EQUATIONS IN D DIMENSIONS 

The field equations of general relativity (in D-dimensional space-time with signature is — I h) form a system of 

ten second-order partial differential equations obeyed by the space-time metric gafi, 

E"P[g,dg,d'g]^^-^T"P[g], (4) 

where the Einstein tensor E"l^ = R"P — ^Rg"P is generated, through the gravitational coupling constant K = StzG/c^, 
by the matter stress-energy tensor T"P . The gravitational constant G is related to the usual three-dimensional Newton's 
constant Gn by 

G = Gw<-^ (5) 

where £o denotes an arbitrary length scale. Among the ten Einstein equations, four govern, via the contracted Bianchi 
identity, the evolution of the matter system, 

^^£"^"=0 =^ v^r"^=o. (6) 

The space-time geometry is constrained by the six remaining equations, which place six independent constraints on 
the ten components of the metric g^p, leaving four of them to be fixed by a choice of a coordinate system. 

In this paper we adopt the conditions of harmonic, or de Donder, coordinates. We define, as a basic variable, the 
gravitational-field amplitude'* 

^^g^P-n'^P, (7) 



^ Here denotes the contravariant metric (satisfying g^^g^p = S^), where g is the determinant of the covariant metric, g = det(gQ^), and where 
rj "'^ represents an auxiliary Minkowskian metric. 



The harmonic-coordinate condition, which accounts exactly for the four equations (|6j corresponding to the conserva- 
tion of the matter tensor, reads 

d^.h"'' = 0. (8) 

The equations Q |S} introduce into the definition of our coordinate system a preferred Minkowskian structure, with 
Minkowski metric rjccp- Of course, this is not contrary to the spirit of general relativity, where there is only one 
physical metric g^ip without any flat prior geometry, because the coordinates are not governed by geometry (so to 
speak), but rather are chosen by researchers when studying physical phenomena and doing experiments. Actually, 
the coordinate condition (|8j is especially useful when we view the gravitational waves as perturbations of space-time 
propagating on the fixed Minkowskian manifold with the background metric rjap ■ This view is perfectly legitimate and 
represents a fruitful and rigorous way to think of the problem when using approximation methods. Indeed, the metric 
rjijip, originaUy introduced in the coordinate condition (|8}, does exist at any finite order of approximation (neglecting 
higher-order terms), and plays in a sense the role of some "prior" flat geometry. 

The Einstein field equations in D dimensions, relaxed by the condition of harmonic coordinates (|8}, can be written 
in the form of inhomogeneous flat d' Alembertian equations, 

ah'^P^^-^T^P, (9) 

where □ = □,)= rj'^^dfidv denotes the D-dimensional flat space-time d' Alembertian operator The source term, t"^, 
can rightly be interpreted as the stress-energy pseudo-tensor (actually, is a Lorentz tensor) of the matter fields, 
described by T"P, and the gravitational field, given by the gravitational source term A"'^, i.e. 

T«/^ = |g|r«/^ + ^A«i^. (10) 

The exact expression of A"'^, including all non-linearities, reads 

- g"^g,rdxhP'd^h'^ ~gP^g,rdxh"'d^h'^ +g^yg''d,h"^d,hP' 

+ \i28''^8^''~8''^8n{8Xr8e.~-^8re8i.)d,h'-''dyh''. (11) 

In this form the only explicit dependence on the dimension D is in the last term of il 1> . As is clear from this expression, 
A"^ is made of terms at least quadratic in the gravitational-field strength h and its first and second space-time 
derivatives. In the following, for the highest post-Newtonian order that we consider (3PN), we need the quadratic, 
cubic and quartic pieces of A"'^ . With obvious notation, we can write them as 

A"P =N"P[h,h]+M"P[h,h,h]+L"P[h,h,h,h] + ff{h^). (12) 



These various terms can be straightforwardly computed from Eq. ( II 1> : see Eqs. (3.8) in Ref. 13511 for explicit 
expressions. 

As said above, the condition is equivalent to the matter equations of motion, in the sense of the conservation of 
the total pseudo-tensor t"P , 

d^T"^ = o ^ v^r«'' = o. (13) 



When developing post-Newtonian approximations, we look for the solutions of the field equations (191 11 OI[TT1ll3t under 
the following four hypotheses: 

1. The matter stress-energy tensor T"P is of spatially compact support, i.e. can be enclosed into some time-like 
world tube, say r <a, where r = |x| is the harmonic-coordinate radial distance. Outside the domain of the source, 
when r > a, the gravitational source term, according to Eq. il3\ . is divergence-free, 

a^A"^=0 (whenr>fl). (14) 

2. The matter distribution inside the source is smooth: T'^P e C°°{W') where li — D — 1 is the space dimension. We 
have in mind a smooth hydrodynamical "fluid" system, without any singularities nor shocks (a priori), that is 
described by some Eulerian equations including high relativistic corrections. 



The source is post-Newtonian in the sense of the existence of the small post-Newtonian parameter v jc— OiXjc). 
For such a source we assume the legitimacy of the method of matched asymptotic expansions for matching the 
inner post-Newtonian field, which is valid only in the source's near zone, and the outer multipolar decomposition 
in the source's exterior near zone. 

The gravitational field has been independent of time (stationary) in some remote past, i.e. before some finite 
instant —S/' m the past, in the sense that 



dt 



whenf<-^. (15) 



The latter condition is a mean to impose (somewhat by brute force), the famous no-incoming radiation condition, 
ensuring that the matter source is isolated from the rest of the Universe and does not receive any radiation from 
infinity. Ideally, the no-incoming radiation condition should be imposed at past null infinity. One can argue f^?] that 
the condition of stationarity in the past ( I15t . although much weaker than the real no-incoming radiation condition, 
does not entail any physical restriction on the validity of the formulas we derive. 

Subject to the condition jl5L the Einstein differential field equations can be written equivalently into the form 
of the integro-differential equations 

/,«/^(x,f) = l^ /t/Vt/f'Gret(x-x',f-f')T"^(x',f')- (16) 



where Gret(x,f) is the scalar retarded Green function in D = c/ + 1 dimensions. The Green function for general d has 
no simple expression in (f,x) space. However, starting from its well-known Fourier-space expression, one can write 
the following simple integral expression (see e.g. |43]), 



0(0 t^-(k 



[i%Yi^ Jo 



GRet(x,f) = -7 — T779 / dk{-] sm{ckt)Jd _Akr). (17) 



Notice that this is in fact a function of t and f" = |x| only. Here 9{t) is the Heaviside step function, and Jd _y{kr) the 
usual Bessel function [see Eq. ( I73t in the Appendix]. 



3. DIMENSIONAL REGULARIZATION OF THE EQUATIONS OF MOTION 

As said before, work at the 3PN order using Hadamard's self-field regularization showed the appearance of ambiguity 
parameters, due to an incompleteness of the Hadamard regularization employed for curing the infinite self-field of 
point particles. By ambiguity parameter, we mean a dimensionless coefficient which cannot be computed within 
the Hadamard regularization scheme. Nevertheless, the majority of terms could be computed unambiguously with 
Hadamard's regularization i34[ Issll . We summarize here the determination using dimensional regularization of the 
ambiguity parameter X which appeared in the 3PN equations of motion; note that X is equivalent to the static ambiguity 
parameter ©static originally introduced in Refs. |44, 45]. 

The post-Newtonian iteration of the Einstein field equations with point-like matter source (Dirac delta-functions 
with spatial supports yi and y2) yields a generic form for the functions representing the metric components in 
successive post-Newtonian approximations. The generic functions we have to deal with in 3 dimensions, say F{x), are 
smooth on M? except at yi and yi, around which they admit singular Laurent-type expansions in powers and inverse 
powers of ri = |x — yi | and r2 = |x — y2|, given by (say, for any A^)^ 

E rf/,(ni)+o(rf), (18) 

po<P<A' 1 

and similarly for the other point 2. Here ri = |x — yi | ^ 0, and the coefficients i/p of the various powers of r\ depend 
on the unit direction ni = (x — y\)/r\ of approach to the singular point. ^ The powers p of r\ are relative integers, and 



' The function f (x) depends also on time (, through for instance its dependence on the velocities vi(?) and V2(f), but the (coordinate) / time is 
purely "spectator" in the regulaiization process, and thus will not be indicated. 
^ The o Landau symbol for remainders takes its standard meaning. 



are bounded from below (pQ < p). The coefficients \ fp (and 2/p) for which p <0 can be referred to as the singular 
coefficients of F. 

A function F defined in 3 dimensions being given, and admitting the singular expansion ( I18t . we define the 
Hadamard partiefinie of F at the location of the particle 1, where it is singular, as 

(ni=/^/o(n,), (19) 

where dQ.\ denotes the solid angle element centered on y 1 and of direction ni . An important feature we have to notice 
is that because of the angular integration in Eq. M9\ , the Hadamard partie finie is "non-distributive" in the sense 
that {FG)\ ^ {F)\{G)\ in general. The non-distributivity of Hadamard's partie finie will be the main source of the 
appearance of ambiguity parameters at the 3PN order; remarkably it does not affect any calculation before that order. 
The second notion of Hadamard partie finie (in short Pf) concerns that of the integral / d^xF, which is generically 
divergent at the location of the two singular points y 1 and y2 (we assume for the moment that the integral converges at 
infinity). It is defined by 

Pf,,,, /t/\F = lim| / d\F + 4n V f— ( ^\ +47rln f-Vr^F), + 1 ^ 2I. (20) 

The first term integrates over a domain S^{s) defined as from which the two spherical balls ri < s and r2 < s of 
radius s and centered on the two singularities are excised. The other terms, where the value of a function at point 
1 takes the meaning ( I19l l. are such that they cancel out the divergent part of the first term in the limit where — > 
(the symbol 1^2 means the same terms but corresponding to the other point 2). The Hadamard partie-finie integral 
depends on two strictly positive constants si and S2, associated with the logarithms present in Eq. ( I20l l. 

The post-Newtonian approximation consists of breaking the hyperbolic d' Alembertian operator □ in Eq. (|9j into 
the elliptic Laplacian operator A plus the famous retardation term —c^'^d} which is to be considered small in the 
post-Newtonian sense, and put in the RHS of the equation where it is iterated. As a consequence, we essentially have 
to deal with the regularization of Poisson integrals, or iterated Poisson integrals, of the generic function F . In the case 
of a Poisson integral potential in 3 dimensions, say 

the Hadamard partie finie integral must be defined in a more precise way. Indeed, the definition ( I19> stricto sensu 
is applicable when the expansion of the function F, when r\ 0, does not involve logarithms of ri\ see Eq. ( I18> . 
However, the Poisson integral P(x') of F(x) will typically involve such logarithms (these will appear precisely at the 
3PN order), namely some Inr'j where r'j = |x' — yi | tends to zero (hence Inr'j is formally infinite). The proper way to 
define the Hadamard partie finie in this case is to include the Inr'j into its definition, and we arrive at I4(tIi 

^F(x)-(r?F)i. (22) 

The first term follows from Hadamard's partie finie integral ilOi : the second one is given by Eq. ( I19I I. Notice that in 
this result the constant si entering the partie finie integral (I20> has been "replaced" by r'j , which plays the role of a new 
regularization constant (together with for the other particle), and which ultimately parametrizes the final Hadamard 
regularized 3PN equations of motion. It was shown that r'j and are unphysical, in the sense that they can be removed 
by a coordinate transformation L34..35 1. On the other hand, the constant ^2 remaining in the result i22\ is the source 
for the appearance of the physical ambiguity parameter, called A, as it will be related to it by Eq. ( I28t below. 

In d spatial dimensions, there is an analogue of the function F, which results from the same post-Newtonian iteration 
process but performed in d spatial dimensions. Let us call this function F^''\x), where x e . When ri ^ the 
function F^''^ admits a singular expansion which is richer than in 3 dimensions. Posing £ = c/ — 3 we have 

f(^)(x)= ^ rr'V|4(ni)+«('f)- (23) 

P0<p<N 1 



(£) 

The coefficients /ip,^(ni) depend on the dimension; the powers of ri involve the relative integers p and q whose 
values are limited by some pQ, qo and qi as indicated. The Poisson integral of F^''\ in d dimensions, is given by the 
Green's function for the Laplace operator, 

P^'\^')--^l'r^^F^'H^l (24) 



(25) 



4n 

where ^ is a constant related to the usual Eulerian F-function by 



d-2 



We have lim^/^g ^ = 1 . Notice also that k is closely linked to the volume Q.^^ \ of the sphere with d~l dimensions; 
see Eq. \b5\ in the Appendix. 

We need to evaluate the Poisson integral at the point x' = yi where it is singular; in contrast with the case of 
Hadamard regularization where the result is given by (I22> . this is quite easy in dimensional regularization, because 
the nice properties of analytic continuation allow simply to get [P'^'^)(x')]x'=y, by replacing x' by yi into the explicit 
integral form MA\ . So we simply have 

P^'Hy,)^~^l^,F^"H-). (26) 

The main technical step of our strategy to compute the ambiguity parameter X will consist of computing, in the limit 
e — > 0, the difference between the t/-dimensional Poisson potential ( I26t . and its 3-dimensional counterpart which is 
defined from Hadamard's self-field regularization and given by i22\ . Denoting the difference between the dimensional 
and Hadamard regularizations by means of the script letter ^, we pose (for the result concerning the point 1) 

^PiEEp('')(yi)~(P)l. (27) 

That is, S^Pi is what we shall have to add to the Hadamard-regularization result in order to get the li-dimensional 
result. However, we shall only compute the first two terms of the Laurent expansion of SPi when e = d — 3 ^ 0, say 
a_i e^' +ao + ^{e). This is the information we need to determine the value of the ambiguity parameter. Notice that the 
difference !^Pi comes exclusively from the contribution of terms developing some poles 1/g in the d-dimensional 
calculation. 

Let us next outline the way we obtain, starting from the computation of the "difference", the 3PN equations of 
motion in dimensional regularization, and show how the ambiguity parameter A is fixe d by the process. By contrast to 
r'j and which are pure gauge, A is a genuine physical ambiguity, introduced in Refs. |46, 35] as the single unknown 
numerical constant parametrizing the ratio between S2 and ^2 [where S2 is the constant left in Eq. J22H as 

ln(^)=i|+A^^ (and 1^2), (28) 
Vs2^ 308 m2 

where mi and m2 are the two masses. The terms corresponding to the A-ambiguity in the acceleration ai — d\\/dt of 
particle 1 read simply 

A^.[;il- 44A G%miml{mi+m2) ^ 

AaiiAJ- -^-g ni2, (29) 

where the relative distance between particles is denoted yi —y2 = '"121112 (with 1112 the unit vector pointing from 
particle 2 to particle 1). We start from the end result of Ref. i35ll for the 3PN harmonic coordinates acceleration ai in 
Hadamard's regularization, abbreviated as HR. Since the result was in fact obtained by means of a specific variant of 
HR, called the extended Hadamard's regularization (in short EHR), we write it as 



(HR) _ (EHR) 



-Aai[A], (30) 



where aj™'^' is a fully determined functional of the masses mi and m2, the relative distance ri2ni2, the coordinate 
velocities vi and V2, and also the gauge constants r'j and rj. The only ambiguous term is the second one and is given 
by Eq. (|29|i. 



Our method is to express both the dimensional and Hadamard regularizations in terms of their common "core" 
part, obtained by applying the so-called "pure-Hadamard-Schwartz" (pHS) regularization. Following the definition 
of Ref. 1 39 1, the pHS regularization is a specific, minimal Hadamard-type regularization of integrals, based on the 
partie finie integral i20\ . together with a minimal treatment of "contact" terms, in which the definition i20\ is applied 
separately to each of the elementary potentials, denoted V, V,-, Wij • • • , that enter the post-Newtonian metric. The 
pHS regularization also assumes the use of standard Schwartz distributional derivatives |47J. The interest of the pHS 
regularization is that the dimensional regularization is equal to it plus the "difference"; see Eq. (I33> below. 

To obtain the pHS -regularized acceleration we need to substract from the EHR result a series of contributions, which 
are specific consequences of the use of EHR L46i. .48J . For instance, one of these contributions corresponds to the fact 
that in the EHR the distributional derivative differs from the usual Schwartz distributional derivative. Hence we define 



(pHS) _„(EHR) 
^1 - ^1 



L^Aai, (31) 

A 

where the S^ai's denote the extra terms following from the EHR prescriptions. The pHS -regularized acceleration j3U 
constitutes essentially the result of the first stage of the calculation of ai, composed of plenty of terms and which are 
all perfectly well-defined. 

The next step consists of evaluating the Laurent expansion, in powers of e = c/ — 3, of the difference between 
the dimensional regularization and the pHS (3-dimensional) computation. As we said above this difference makes a 
contribution only when a term generates a pole ~ 1/e, in which case the dimensional regularization adds an extra 
contribution, made of the pole and the finite part associated with the pole [we consistently neglect all terms &"{£)]■ 
One must then be especially wary of combinations of terms whose pole parts finally cancel ("cancelled poles") but 
whose dimensionally regularized finite parts generally do not, and must be evaluated with care. We denote the above 
defined difference by 

2iai=Y^2iPi. (32) 

It is made of the sum of all the individual differences of Poisson or Poisson-like integrals as computed in Eq. J27> . The 
total difference ( 13 2> depends on the Hadamard regularization scales r'j and S2 (or equivalently on A and r'j, rj), and on 
the parameters associated with dimensional regularization, namely e and the characteristic length scale Iq introduced 
in Eq. Finally, our main result is the explicit computation of the dimensional regularization (DR) acceleration as 

ar*=a(P"') + ^ai. (33) 
With this result we can prove two theorems ll39il . 

Theorem 1 The pole part 1/e of the DR acceleration \33\ can be re-absorbed (i.e., renonnalized) into some shifts 
of the two "bare " world-lines: yi —>■ yi + and y2 ^ y2 + ^ith, say, 4 i 2 °^ l/£, so that the result, expressed in 
terms of the "dressed" quantities, is finite when £ — !■ 0. 

The precise shifts and ^2 involve not only a pole contribution 0= 1/e [which would define a renormalization by 
minimal subtraction (MS)], but also a finite contribution when e ^ 0. Their explicit expressions read: 
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1 21n^'^^^'M 327 



e \ ^0 / 1540 



Sni (together with 1 ^ 2), (34) 



where Gn is Newton's constant, £0 is the characteristic length scale of dimensional regularization [cf. Eq. Q], aA^i is 
the Newtonian acceleration of the particle 1 ind dimensions, and q = 4-7te^ depends on Euler's constant C = 0.577 • • • . 
Note that when working at the level of the equations of motion (not considering the metric outside the world-lines), 
the effect of shifts can be seen as being induced by a coordinate transformation of the bulk metric as in Ref. Issll . 

Theorem 2 The renormalized (finite) DR acceleration is physically equivalent to the Hadamard- regularized (HR) 
acceleration (end result of Ref. k3^1 ). in the sense that 

(HR) ,. r (DR) , 5. 1 

a^j ' = lim[ai ' + 5^aiJ, (35) 

where 5^ a\ denotes the effect of the shifts on the acceleration, if and only if the HR ambiguity parameter X entering 
the harmonic-coordinates equations of motion takes the unique value 



The coefficient A has therefore been fixed by dimensional regularization, both within the ADM-Hamihonian formal- 
ism 1 38], and the harmonic-coordinates equations of motion |39]. 

An alternative work, by Itoh and Futamase |49, 50], following previous investigations in Refs. ll5lll52il . has derived 
the 3PN equations of motion in harmonic coordinates by means of a variant of the famous "surface-integral" method 
introduced by Einstein, Infeld and Hoffmann ] 1 ] . The aim is to describe extended relativistic compact binary systems in 
the strong-field point particle limit ] 53 ] . This approach, alternative to the use of self-field regularizations, is interesting 
because it is based on the physical notion of extended compact bodies in general relativity, and is free of the problems 
of ambiguiti^. The end result of Refs. ll49l Isoll is in agreement with the 3PN harmonic coordinates equations of 
motion Il34li35il and, moreover, it does determine the ambiguity parameter X to exactly the value J36> . 



4. DIMENSIONAL REGULARIZATION OF THE RADIATION FIELD 

We now address the similar problem concerning the binary's gravitational radiation field (3PN beyond the Einstein 
quadrupole formalism), for which three ambiguity parameters, denoted k and ^, have been shown to appear due to 
the Hadamard self-field regularization P6, 37]. To apply dimensional regularization, we must use as we did for the 
equations of motion in Section |3] the ^/-dimensional post-Newtonian iteration; and, crucially, we have to generalize 
to d dimensions some key results from the gravitational wave generation formalism. The specific wave generation 
formalism we employ is based on a post-Newtonian expansion for the metric field in the near zone of the source, 
and on the so-called multipolar-post-Minkowskian expansion for the field in the exterior of the source, including 
the regions at infinity from the source where the detector is located ]54]. The expression of the multipole moments 
describing the physical (post-Newtonian) source are then obtained by a technique of asymptotic matching performed 
in the overlapping region of common validity between the two types of expansion, namely the exterior part of the near 
zone L5S. 561 (see ]42] for a review). 

Let us first recall the expressions of the source multipole moments II (mass-type moment) and Jl (current-type) of 
an isolated post-Newtonian source in ordinary 3-dimensional space. They are given, for multipolarities £ > 2, by \56\' 



2(2^+1) 1 
^ c\l+m + 2){U + 5) ^^'+2 (x,f +z|x|/c), (37a) 

2i+l .1 

Se+i{z)xL_i>ac'^hc Mx,f + z|x|/c), (37b) 



c^{£ + 2){2£ + 3) 

where the source densities E^y's are evaluated at the position x and at time f +z|x|/c, and are defined by 

c-^ c 

Here t"I^ is the pseudo stress-energy tensor ilOi in 3 dimensions and the overbar refers to its formal post-Newtonian 
expansion, t"^ = PN[t"'^]. Let us note that the expressions J37> are "exact", in the sense that they are formally valid 



' Our notation is tlie following: L = r'l 12- ■■ ii denotes a multi-index, made of / spatial indices; similarly we write for instance aL = ai[ ...ii or 
L — 1 = I'l . . . . The symmetric-trace-free (STF) projection is denoted with a hat, so that = x^^^ is the STF projection of the product of / spatial 
vectors, denoted xt = x, j . . . .v;, . Sometimes we also indicate the STF projection by brackets surrounding indices xt = xi^u^ . Note that an expansion 
into STF tensors = xi^jr^ (which are functions of the spherical angles 9 and (^) is equivalent to the usual expansion in spherical harmonics 
Y;m(Si0)- The dots indicate successive time-derivations. 



up to any PN order. Equations (I37> involve an integration over the variable z, with associated function 5({z) given by 



5e{z) 



(2^+1)! 
2^+1^1 



2\t 



dzdfXz) 



lim 5i{z) = 



(39) 



where 5{z) is the usual Dirac's one-dimensional delta-function. In practice, the post-Newtonian moments i31i are to 
be computed by means of the infinite post-Newtonian series 



dz 5i{z)T.ix,t+z\x\/c) ■ 



1^ (2^+1)!! 



1)! 



\x[d_ 
c dt 



2k 



E(x,f) 



(40) 



In Eqs. (I37t there is a special process of taking the "finite part", indicated by the symbol J?"^^, which is necessary 
in order to deal with the bound of the integral at infinity, corresponding to infra-red divergencies, in the limit |x| °o. 
Indeed, notice that the pseudo stress-energy tensor t"I^ includes the crucial contribution of the gravitational field, 
denoted A"/^ in Eq. ( II Ot . which has a spatially non-compact support. This fact, together with the presence of the 
multipolar factor xl in the integrand, prevents one to (naively) write the standard expressions for the multipole 
moments valid for compact-supported sources; such expressions have no meaning in non-linear general relativity. 
The solution to this dilemna has been to introduce I5xl5al the specific finite part and to show how these specific 
multipole moments so defined are related to the physical wave form at infinity. 

To proceed with dimensional regularization, we need the cZ-dimensional analogues of the multipole moments of 
the source, say if andjf , consequences of the D-dimensional Einstein field equations for isolated post-Newtonian 



sources. In the case of the mass-type moments we find |4] 



2{d-2) J I m c^{d + l-2){d + 2l) ""-[i+iP ' ' 

2(^ + 2/ -2) „ ^ 

^c\d + l~l){d + l-2){d + 2l + 2) -""'^ [Ctir""' ^> 



where we denote [generalizing Eq. j38H 

^_ 2 (^-2)tQ° + t'' 
d^ c^ ' 



E t'^' 



(41) 



(42) 



where now x"^ is the post-Newtonian expansion of the pseudo stress-energy tensor in D dimensions. For any of the 
latter source densities the underscript [I] means the infinite series 



+^ 1 r(f + /) 



(43) 



which constitutes the cZ-dimensional version of Eq. i40i . At Newtonian order Eq. ( 14 H reduces to the standard result 
l['') = Jd''xpxL + with p = T'^/c^. 

Like for the case of the equations of motion, the ambiguity parameters £, , K and ^ come from a deficiency of the 
Hadamard regularization coming up at 3PN order and mainly due to its "non-distributivity". The terms corresponding 
to these ambiguities are contained in the 3PN mass quadrupole moment 1,-^ and are found to be 



AI,v[^,'c,C] 
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(44) 



where yi, vi and ai denote the first particle's position, velocity and acceleration (the brackets () surrounding indices 
refer to the STF projection). As in Section|3] we express the Hadamard and dimensional results in terms of the more 
basic pure-Hadamard-Schwartz (pHS) regularization. The first step of the calculation L37ijl is therefore to relate the 



HR) 



Hadamard-regularized quadrupole moment !■■ , for general orbits, to its pHS part; we find: 
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(45) 



In the RHS we see both the pHS part, and the effect of adding the ambiguities, with some numerical shifts of the 
ambiguity parameters coming from the difference between the specific Hadamard-type regularization scheme used in 
Ref. 1 36] and the pHS one. The pHS part is free of ambiguities but depends on the gauge constants r'j and introduced 
in the harmonic-coordinates equations of motion |34, 35]. 

We next use the J-dimensional moment J41l i to compute the difference between the dimensional regularization 
(DR) result and the pHS one ] 14, 40]. As in the work on equations of motion, we find that the ambiguities arise solely 
from the terms in the integration regions near the particles {i.e., ri = |x — yi | ^ or r2 = |x — y2| ^ 0) that give 
rise to poles °= 1/e, corresponding to logarithmic ultra-violet divergences in 3 dimensions. The infra-red region at 
infinity {i.e., |x| +°o) does not contribute to the difference DR — pHS. The compact-support terms in the integrand 
of (14 1> . proportional to the components of the matter stress-energy tensor T"P, are also found not to contribute to the 
difference. We are therefore left with evaluating the difference linked with the computation of the non-compact terms 
in the expansion of the integrand in (I41> near the singularities that produce poles in d dimensions. 

Let F^''\x) be the non-compact part of the integrand of the quadrupole moment J41t (with indices L = ij), where 
f''^^ includes the appropriate multipolar factors such as so that 



I]"' = J d^xF^^^x). (46) 

We do not indicate that we are considering here only the non-compact part of the moments. Near the singularities 
the function F^''' (x) admits a singular expansion of the type ( I23> . In practice, the coefficients ifj,^^ are computed by 
performing explicitly the post-Newtonian iteration. On the other hand, the analogue of Eq. ( I46> in 3 dimensions is 

l,j = Pf J d\F{x), (47) 

where Pf refers to the Hadamard partie finie defined by Eq. ( 12m . The difference f^I between the DR evaluation of the 
c/-dimensional integral ( I46L and its corresponding three-dimensional evaluation, i.e. the partie finie ( 147 L reads then 

(48) 

Such difference depends only on the ultra-violet behavior of the integrands, and can therefore be computed "locally", 
i.e. in the vicinity of the particles, when ri — > and — > 0. We find that Eq. (I48t depends on two constant scales si and 
S2 coming from Hadamard's partie finie ( I20l i. and on the constants belonging to dimensional regularization, which are 
e = d — 3 and the length scale £o defined by Eq. (|5}- The dimensional regularization of the 3PN quadrupole moment 
is then obtained as the sum of the pHS part, and of the difference computed according to Eq. ( I48> . namely 

j(DR)^j(pHS)^^j..^ (49) 

An important fact, hidden in our too-compact notation i49i . is that the RHS of ( I49> does not depend on the Hadamard 
regularization scales si and S2, which cancel out from the two terms in the RHS. Therefore it is possible to re-express 
these two terms (separately) by means of the constants r'j and ri, instead of si and S2, where r'j, r2 are the two fiducial 
scales entering the Hadamard-regularization result (I45> . This replacement being made the pHS term in Eq. ( I49> is 
exactly the same as the one in Eq. ( B3t . At this stage all elements are in place to prove the following theorem 



Theorem 3 The DR quadrupole moment H491 is physically equivalent to the Hadamard-regularized one ( end result of 
Refs. I 36n37il ). in the sense that 

(50) 



l(f) = lim 



T DR) c. , 



where 5^\jj denotes the effect of the same shifts as determined in Theorems^^\ if and only if the HR ambiguity 
parameters B,, K and C, take the unique values 

9871 7 



Moreover, the poles 1/e separately present in the two terms in the brackets of \5(M cancel out, so that the physical 
(renormalized or "dressed") DR quadrupole moment is finite and given by the limit when e ^0 as shown in Eq. i l56>t . 



This theorem finally provides an unambiguous determination of the 3PN radiation field by dimensional regularization. 

It should be emphasized that though the values i5H represent the end result of dimensional regularization, several 
alternative calculations have provided a check, independent of dimensional regularization, for all the parameters i5H . 
hi ll37il we computed the 3PN binary's mass dipole moment I,- using Hadamard's regularization, and identified I, 
with the 3PN center of mass vector position G;, already known as a conserved integral associated with the Poincare 
invariance of the 3PN equations of motion in harmonic coordinates |57]. This yields ^ + K = —9871/9240 in 
agreement with Eq. ( I51> . Next, we considered L58.1 the limiting physical situation where the mass of one of the 
particles is exactly zero (say, OT2 = 0), and the other particle moves with uniform velocity. Technically, the 3PN 
quadrupole moment of a boosted Schwarzschild black hole is computed and compared with the result for 1,^ in the 
limit m2 = 0. The result is = —7 /33, and represents a direct verification of the global Poincare invariance of the wave 
generation formaUsm. Finally, fc = is proven i40ll by showing that there are no dangerously divergent "diagrams" 
corresponding to non-zero fc-values, where a diagram is meant here in the sense of Ref. |59]. All these verifications 
confirm the validity of dimensional regularization for describing the dynamics of systems of compact bodies. 



5. CONCLUSION 

The determination of the values (I36I I and ( 15 11 1 completes the problem of the general relativistic prediction for the 
templates of inspiralling compact binaries up to 3PN order, and actually up to 3.5PN order as the corresponding "tail 
terms" composing this order have already been determined ll60tl . The relevant combination of the parameters \5\l 
entering the 3PN energy flux in the case of circular orbits, namely 9 f3^, is now fixed to 

1 1 1 

= <^+2)c+C = --— — . (52) 
^ ^ 9240 

The orbital phase of compact binaries, in the adiabatic inspiral regime {i.e., evolving by radiation reaction), involves at 
3PN order a linear combination of Q and of the equation-of-motion related parameter X 1 13], which is determined as 

7 , 1039 

= - -A = . (53) 

3 4620 

The parameter A appears here because the orbital phase follows from energy balance between the total radiated energy 
flux and the decrease of orbital center-of-mass energy which is computed from the equations of motion. 

The practical implementation of the theoretical templates in the data analysis of detectors follows the standard 
matched filtering technique. The raw output of the detector o{t) consists of the superposition of the real gravitational 
wave signal /ireai(0 and of noise n{t). The noise is assumed to be a stationary Gaussian random variable, with zero ex- 
pectation value, and with (supposedly known) frequency-dependent power spectral density 5„(a)). The experimenters 
construct the correlation between o(f ) and a filter q{t), i.e. 

c{t) = p dt'o{t')q{t + t'), (54) 

and divide c(f) by the square root of its variance, or correlation noise. The expectation value of this ratio defines 
the filtered signal-to-noise ratio (SNR). Looking for the useful signal /!reai(0 in the detector's output o(f), the 
experimenters adopt for the filter 

m = 1^, (55) 

S„{(0) 

where q{(i)) and h{(x)) are the Fourier transforms of q{t) and of the theoretically computed template h{t). By the 
matched filtering theorem, the filter \55\ maximizes the SNR if hit) = /!reai(0- The maximum SNR is then the best 
achievable with a linear filter In practice, because of systematic errors in the theoretical modelling, the template h{t) 
will not exactly match the real signal /ireai(0' but if the template is to constitute a realistic representation of nature 
the errors will be small. This is of course the motivation for computing high order post-Newtonian templates, in 
order to reduce as much as possible the systematic errors due to the unknown post-Newtonian remainder. The fact 
that the numerical value of the parameter \53\ is quite small, Q ~ 0.22489, indicates, following measurement-accuracy 
analyses [lOi 111, that the 3PN (or better 3. 5PN) templates should constitute an excellent approximation for the analysis 
of gravitational wave signals from inspiralling compact binaries. 



ACKNOWLEDGMENTS 



It is a pleasure to thank Thibault Damour, Gilles Esposito-Farese and Bala Iyer for the exciting collaborations f3 
summarized here, and to warmly acknowledge Jean-Michel Alimi for organizing a very interesting conference. 



A. USEFUL FORMULAS IN d SPATIAL DIMENSIONS 

This Appendix is intended to provide a compendium of (mostly well-known) formulas for working in a space with 
d dimensions. As usual, though we shall motivate some formulas below by writing some intermediate expressions 
which make complete sense only when is a strictly positive integer, our final formulas are to be interpreted, by 
complex analytic continuation, for a general complex dimension, e C. Actually one of the main sources of the 
power of dimensional regularization |29, 30!. l3ll 13211 is its ability to prove many results by invoking complex analytic 
continuation in d. 

We discuss first the volume of the sphere having d — \ dimensions, i.e. embedded into Euclidean li-dimensional 
space. We separate out the infinitesimal volume element in d dimensions into radial and angular parts, 

d''x^r''-^drdQ.j_u (56) 



where r = |x | denotes the radial variable (/. e. , the Euclidean norm of x e K ) and dQ.d- 1 is the infinitesimal solid angle 
sustained by the unit sphere with d — I dimensional surface. To compute the volume of the sphere, = / dQ.d-\, 
one notices that the following li-dimensional integral can be computed both in Cartesian coordinates, where it reduces 
simply to a Gaussian integral, and also, using i56\ . in spherical coordinates: 

d'xe-'-' ^l^j dxe-'^'y' ^ ni ^ dr^^'-'e'^ = ^a,_iT (^0 , (57) 

where F in the last equation denotes the Eulerian function. This leads to the well known result 

For instance one recovers the standard results Q.2 = 4n and £2i = 2n, but also Ho = 2, which can be interpreted by 
remarking that the sphere with dimension is actually made of two points. If we parametrize the sphere i in li — 1 
dimensions by means of li — 1 spherical coordinates O^-i, Qd-ii ' " ^ which are such that the sphere in d — 2 

dimensions is then parametrized by 9ci-2^ • • s and so on for the lower-dimensional spheres, then we find that the 
differential volume elements on each of the successive spheres obey the recursive relation 

da,,-i - {smed-iY'^dOd-ida,,-!- (59) 



^d-i r.. ^d-2_ 



Note that this implies 

"1' ' = I ded-i{sme,,_i)''-^= I dxil-x"-)^, (60) 

Jo J-l ' 

which can also be checked directly by using the explicit expression ( 15 8> . 

Next we consider the Dirac delta-function 5''''(x) in d dimensions, which is formally defined, as in ordinary 
distribution theory [47.1 . by the following linear form acting on the set of smooth functions G C°°(M'') with compact 
support: yep e ^, 

<d^''\(p>= [ d''xd^''\x)(p{x)^(p{0), (61) 



where the brackets refer to the action of a distribution on (p E S^. Let us now check that the function defined by 

M = kr^-'', (62a) 
k = (62b) 



where r is the radial coordinate in d dimensions, such that — Lf=i(-^')^' Green's function of the Poisson 

operator, namely that it obeys the distributional equation 

Au = -4-nd^''\x). (63) 

For any a e C we have Ar" = a{a + d — 2) r"^^, thus we see that Am = in the sense of functions. Let us formally 
compute its value in the sense of distributions in x-space.** We apply the distribution Am on some test function (jO G ^, 
use the definition of the distributional derivative to shift the Laplace operator from u to (p, compute the value of the 
(/-dimensional integral by removing a ball of small radius s surrounding the origin, say B{s), apply the fact that Am = 
in the exterior of B{s), use the Gauss theorem to transform the result into a surface integral, and finally compute that 
integral by inserting the Taylor expansion of (p around the origin. The proof of Eq. ( l63l is thus summarized in the 
following steps: 

<Am,^> = <M,A^> 

= lim / d''xuA(p 

~ lim / d'^^xdj[udj(p — diU(p] 

— lim / s''^^dQ.d-i{—ni)[udi(p — diU(p] 

= lim / s''-UQ.d-ii~ni) \-~k{2-d)s^-''ni(p{0) 

s^O J I 

= -47r(p(0). (64) 
In the last step we used the relation between k and the volume of the sphere, which is 

An 

kQd-i = -1—,- (65) 
d — 2 

From M = kr^^'^ one can next find the solution v satisfying the equation Av = m (in a distributional sense), namely 



2{4-d)' 



(66) 



From ( I66l l we can then define a whole "hierarchy" of higher-order functions w, ■ ■ ■ satisfying the Poisson equations 
Aw — V, ■ ■ ■ in a distributional sense. 

However, the latter hierarchy of functions m, v, • • • is better displayed using some different, more systematic notation. 

This leads to the famous Riesz kernels, here denoted 5a\ in (/-dimensional Euclidean space j33|-^ These kernels 
depend on a complex parameter a G C, and are defined by 

^''^(x) = Kar"-'', (67a) 

r(^) 

2«;r^r(f) 

For any a G C, and also for any d eC, the Riesz kernels satisfy the recursive relations 

Ad^:%^-d^J\ (68) 



The usual verification of 1631 is done in Fourier space. 
' Besides the Euclidean kernels Sa \ we also have the Minkowski kernels (denoted z'f^) I33II . which are at the basis of the Riesz analytic 
continuation method (2^. 



Furthermore, they obey also an interesting convolution relation, which reads simply, with the chosen normalization of 
the coefficients Ka, as 

(69) 

When a = we recover the Dirac distribution in d dimensions, 8^^^ = K^r^^' = (the coefficient vanishes in this 

case, Kq — 0), and we have u ~ 47152^''', v = —47z5^^\ ■ ■ ■ . 

The convolution relation ( I69> is nothing but an elegant formulation of the Riesz formula in d dimensions. To check 
it let us consider the Fourier transform of r" in d dimensions. 



fa{k)= I d'xlxl'^e-''^ (70) 

Using <56l l we can rewrite it as 



/„{k)=y^ drr"+^-^ J d£ld-ie-"'^ (71) 

in which the angular integration can be performed as an application of Eq. i59\ . This yields an expression depending 
on the usual Bessel function. 



daz-ie-"^" = ^d-i idOd-x (sin0,,_i)''"'e"''''™''''-' - (27r)^ {kr)'-^- J ,^_^{kr), (72) 



where k = |k|, and where we adopt for the Bessel function the defining expression 



r(v + i)r(i) 



■^v W - ' n^/n / (1 ~n (73) 



The radial integration in Eq. J71> is then readily done from using the previous expression, and we obtain 

/„(k)=2«+'';ri-)^^-«-'', (74) 

where the factor in front of the power k^"^'', say A^, is checked from the Parseval theorem for the inverse Fourier 
transform, which necessitates that AaA_a_^/ — {iTt)'' . To obtain Eq. ( 1741 we employ the integration formula 



r 



l+M+v 

dzz^Jv (z) - 2'' ), ^ ( . (75) 

r(i^ 

Finally we can check the Riesz formula by going to the Fourier domain, using the previous relations. The result, 

^-^^-^ r(-f)r(-f)r(H±i±^) ' '''' 



is equivalent to Eq. (I69> . 
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